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Abstract 



Using the Density Matrix Renormalization Group method, we determined 
the phase diagram of a frustrated antiferromagnetic spin-^ ladder at zero 
temperature. Two spin-gapped phases, the Haldane phase and the singlet 
phase, are identified. A phase transition between the two phases occurs at 
any non-zero value of frustration coupling J x . On the phase boundary, the 
spin gap vanishes for sufficiently small J x • Crossing this non-gapped line, 
the transition is of second order, while of first order for larger J x . Striking 
frustration effects are predicted for ladder materials. 

PACS. 75.10.Jm, 75.30.Kz, 75.45.+j 
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In recent years, a lot of theoretical and experimental effort has been devoted to the 
understanding of spin-| antiferromagnetic ladders. fl|-|9| These systems, which interpolate 
structures between one and two dimensions, are realized in materials such as Sr n _ 1 Cu n+ i0 2 n 
HJ. The ladders, being arrays of coupled chains, can display significantly different physical 
properties according to their number of legs. Especially, the spin correlation function in the 
ground state shows an algebraic decay if the number of legs is odd and an exponential decay 
if it is even. j|,[7| This property essentially depends on whether or not the system exhibits 
a spin gap in the spectrum. Therefore, it seems that the ladders are analogous to spin-S* 
antiferromagnetic Heisenberg chains where an odd half-integer spin corresponds to an odd 
number of legs and an integer spin to an even number of legs. [|I0| -|T5|| 



The aim of this paper is to investigate the following two issues: i) How far is the above 
analogy valid? In particular, are the low-lying energy properties of the two-leg standard 
ladder identical to those of the single S = 1 chain? It is clear that a two-leg ladder with a 
ferromagnetic interchain coupling Jj_ < is equivalent to the S = 1 chain when Jj_ — > — oo. 
In the case of the antiferromagnetic interchain coupling, there is also a gap opening at 



any J±, however no phase transition occurs in the whole parameter range. [[RJ. For large 
J_i_, the gap stems from the formation of a spin singlet on each rung. Therefore the origin 
of the spin gap for Jj_ > seems distinguishable from that in the S — 1 chain. [|TJ ii) 
What are the frustration effects caused by a next-nearest-neighbour interaction J x ? Under 
some circumstances, frustration can produce either a dimerized state or a non-dimerized 
spin liquid state. |L7| In the ladder material SrCu203, the frustration is presumed to be 
negligible and J± comparable to the intrachain coupling Jm. |l| In this case, compared to 
the theoretical prediction of A = 0.5J|| [@,|J, the gap is estimated 35% smaller by fitting the 
spin susceptibility and 4% larger from the nuclear magnetic resonance measurement. The 
origin of this discrepancy between the two estimates is unclear. |TJ Our calculation will show 
that for J± = J|| the spin gap only slightly changes as < J x < 0.4J||. The maximum value 
of A is 0.519J|| at J x = 0.25J||. Therefore, it is an open question whether ladder materials 
are sufficiently well described by the standard spin ladder. 

In this work, we study a two- leg frustrated antiferromagnetic spin-| ladder described by: 

H = }^ [Ji^S^j • Si^+i + S2,i • S2,i+i) + J±Si t i ■ S>2,i 

i=l,N 

+ Jx(Sl,i ' S 2 ,i+1 + Si t i+1 ' ^2,*)], (1) 

where S nj j denotes a spin-| operator at site i of the nth chain. J\\ is an intrachain coupling 
between two neighboring spins in each chain, J± an interchain coupling between two spins 
on each rung and J x an interchain coupling between two spins of neighboring rungs. Since 
J x = J\\ is a symmetric line in the parameter space, we only consider the case of J x < Jm. 
When J x = 0, H reduces to that for the standard spin ladder. Hereafter we set J\\ = 1. 

To overcome the inherent difficulties related to the frustration, which plague the Monte 
Carlo simulation, we employ the density matrix renormalization group method invented by 
White. [0 This method has proven to be very powerful and efficient for a systematic study 
of low-lying energy properties of low-dimensional lattice models. In our calculation, we keep 
at least 240 states, and up to 450 states are necessary for critical regimes. The truncation 
error is typically of the order of 10~ 8 . Lengths up to 550 rungs are considered for open 
boundary conditions and finite size scaling is used to determine the thermodynamic limit. 
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For convenience, the number of rungs N is chosen to be even. To check the accuracy and 
the convergence, we performed calculations keeping up to 600 states. The errors on physical 
quantities are estimated to be less than one percent in most cases. The technical details will 



be given elsewhere. [19 



In Fig. 1, we show the phase diagram characterized by two fixed points: (i) the "Haldane 
phase" so named as it contains the limiting case J x = 1 and J± = 0, 

H = Y,(Ki + S 2 ,i) • (Si,+i + S 2 , i+ i) (2) 

whose low-lying spectrum is identical to that of a S = 1 chain; (ii) the "singlet phase" as 
it contains the case J± ^> 1, with a ground state consisting of a singlet on each rung and 
low-lying excitations generated by creating triplets on rungs. In both cases, the system is 
gapped. However, a phase transition (shown in Figs. 3-5) occurs in the parameter space as 
we cross from (i) to (ii). 

Of special interest in the characterization of the two phases, is the sensitivity of their 
low-lying energy spectra to boundary conditions or impurity effects. In particular, both 
experiment ]20|] and numerical calculations [pT] , p2[ indicate that boundary or impurity effects 



are a characteristic feature of the S = 1 chain . 

In the Haldane phase, for open boundary conditions, the ground state is four-fold degen- 
erate with total spin St — 0,1 and the continuum starts with St = 2 as shown schematically 
in Fig. 2a D , while for periodic boundary conditions, the ground state becomes a singlet and 
the lowest excitation has S t — 1 (Fig. 2a p ). We can understand this difference by consider- 
ing the open boundary conditions as a special case of a bond impurity in the weak or strong 
coupling limit p2| , p4| . 

Furthermore, for open boundary conditions only, the region C% appears inside the Hal- 
dane phase, characterized by the presence of midgap states. The number of levels is one as 
we enter the C\ region from below (Fig. 2b) and grows larger and larger as we approach the 
phase boundary (Fig. 2c). This effect is a peculiar feature of this frustrated spin ladder in 
contrast to the S = 1 chain. 

In the singlet phase, independently of the boundary conditions, the ground state is a 
singlet and the continuum starts from a state with S t = 1 outside the region C*2 and S t = 
inside the region C*2 (Fig. 2d-f). In contrast to the Haldane phase, neither midgap states 
for open boundary conditions nor any impurity state due to a bond impurity are found in 
this phase. Q 

On the phase boundary, there are at least three non-trivial degenerate states. Two of 
them have St — and one has S t — 1. As J x < A = 0.287, a non- gapped line shows 
up on the boundary curve (Fig. 2g). Below we will discuss the phase diagram as it is 
numerically determined by the calculation of the ground state energy, singlet density per 
rung and low-lying excitations. 

Fig. 3 shows the ground state energy per rung e . For each J x , a phase transition 
between the Haldane phase (left) and the singlet phase (right) takes place at the maximum 
value of eo- When J x is sufficiently large, eo is obviously singular at the critical value of J±. 
When J x = 1, the singularity originates from the crossing of the two lowest energy levels 
with S t = 0. However, such a crossing disappears due to level repulsion for J x < 1. The 
two lowest energy levels with S t = become degenerate at the critical values only in the 
thermodynamic limit. When J x is small, the curvature of eo looks smooth at its maximum 
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value. The critical values J± t o of J± are given in Table I, which are splined into a solid curve 
to indicate the phase boundary shown in Fig. 1. In the construction of the phase diagram, 
we also verified the transition by varying J x for given Jj_. The results for J± = 1 is shown 
in inset. 

In order to further investigate the ground state properties, we calculate the singlet density 
per rung p s . For a state \^f), we define p s = ^ J2i(^\<SiSj\^) , where Si = ^(| Tii2>i — I I1T2 
)i) is the singlet state formed by two spins on the ith rung. For J x = 1 and in low-lying 
energy states, p s = 1 when Jj_ > Jj_ j0 ; otherwise p s = 0. For J x < 1, we show p s for the 
ground state in Fig. 4. At J± t o, p s changes abruptly for large J x and smoothly for small 
J x . We note that p s = \ at J± = J x = 0, p s < 7 in the Haldane phase, and p s > \ in the 
singlet phase. Therefore, p s also characterizes all the phases. 

We now discuss the low-lying excitations which are used to determine the regions C\ 
and C2 as well as to provide evidence for the existence of a non-gapped line. The spin 
gap A governs the asymptotic behavior of the correlation function in the ground state and 
the low-temperature behavior of the thermodynamic quantities. A is defined as an energy 
difference between the ground state and the bottom state of continuum for spin systems. In 
Fig. 5, we show the energies of several relevant states relative to the ground state energy. 

For J x = 1, the total spin on each rung is conserved. Consequently, the states with n- 
singlets on rungs are the same as those of a S = 1 chain doped by n non-magnetic impurities 
at sites corresponding to those rungs at which the n-singlets are localized. In the Haldane 
phase, the midgap states are the lowest impurity states for the given number of singlets or 
impurities. For instance, two lowest midgap states shown in inset are the ground states of 
open S = 1 chains with N — 1 and N — 2 spins: the lowest one involves a singlet at one 
end of an open ladder and the other two singlets at both ends. In the singlet phase, the 
low-lying excitations are discrete but infinitely degenerate. Two lowest excitations shown in 
inset are generated from the ground state by introducing one, two, and four neighbouring 
triplets on rungs. Their St equal to one when one triplet is involved, or else St = 0. 

For J x < 1, since the total spin on each rung is no longer conserved, the low- lying energy 
properties are dramatically changed according to the magnitude of J x and J±. As shown 
in Fig. 5, there is no midgap state at J x = 0.2 and 0.4. For a larger J x , midgap states 
show up when J± is larger than a critical value Jx,i- Splining J x and J± t i given in Table I, 
we obtained a critical values for J x of 0.430. In the singlet phase, the infinite degeneracy 
at J x = 1 is removed by infinitesimal J11 — J x so that the discrete levels become different 
bands. The bottom state has S t equal to 1 or depnding on J x and J±. For sufficiently 
large J x , the two lowest excitations with S t = 0, 1 across at J± = J|_ i2 . Several J_|_ >2 given in 
Table I are used to determine the value of J x , namely A = 0.287, at which the splined curve 
meets the phase boundary. We note that this leads to a critical point in the phase diagram 
because the splined curve has a curvature clearly opposite to that of the phase boundary. 

As seen from Fig. 5, for each J x < 1 the spin gap has a minimum given in Table I at 
Jfo- For sufficiently large J x , A has a jump at Jj_ j0 , e.g. see Fig. 5 for J x = 0.8. We 
found that the jump and the minimum value, and the width of C 2 decrease coherently as J x 
becomes small and approach to zero for J x — > A. Therefore, a non-gapped line exists on the 
phase boundary for < J x < A. On this line, p s = 1/4 and the ground state has a different 
symmetry than in the other phases. |B1 Further taking into account the critical behavior of 



e and p s , we conclude that the transition between the Haldane phase and the singlet phase 
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crossing the non-gapped line is of second order, while of first order for J x > A. 

In Fig. 6, we show the behavior of the spin gap varying J x for Jj_ = Jii = 1, as is relevant 
to ladder materials such as SrCu203, It is remarkable that the spin gap is almost constant 
when < J x < 0.4. At J x = 0.25, it has a maximum value of A = 0.519. Accordingly, it 
would be interesting to examine the relevance of frustration in ladder materials. 

We have established the phase diagram of a frustrated spin-| ladder using the low-lying 
energy properties. Clearly, the low-lying energy properties of the two-leg antiferromagnetic 
spin- 1 ladder is different from those of the S = 1 chain. According to our findings, one can 
naturally expect similar effects for a frustrated t — J ladder. In addition, for very low carrier 
doping, bound states should occur in the Haldane gap. As far as materials are concerned, 
it would be interesting to explore two kinds of impurity effects: non-magnetic and magnetic 
doping, which produce midgap states for sufficiently large value of J x and impurity states 
for any non-zero value of J x in the Haldane gap, respectively. Moreover, in order to observe 
the phase transition, one could measure the specific heat and spin susceptibility on pure 
ladder materials for different interchain coupling. 

The author is very grateful to Martin Long, Felix Naef and Xenophon Zotos for their 
kind help and stimulating discussions. He also thanks A. Honecker, I. Peschel, and Y. Zhao 
for discussions. This work is supported by the Swiss National Fond Grant No. 20-49486.96. 
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TABLES 



TABLE I. For several values of J x , J±,o, J±,i, J±,2 are the values of J± for points on the phase 
boundary, lower boundary of C\ and upper boundary of C2, respectively; A m j n is the minimum of 
A. Lower numbers are errors at the last digit. 



Jx 


0.2 


0.35 


0.4 


0.5 


0.6 


0.8 


1 


J±,o 


0.378 3 




0.715 2 




I.OOO5 


1.225 5 


1.4015 


J±,l 








0.792 8 


0.835 5 


0.917 5 


0.991 


J±,2 




0.650s 


0.742 8 


0.921s 


LH83 


1.560 5 


2 




0.004 4 




0.027 5 




0.1305 


0.3182 


0.4105 
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FIGURES 

FIG. 1. Phase diagram: a thick solid curve splining dots indicates the phase boundary; Dotted 

line the symmetric line; C\ is surrounded by the phase boundary, a thin curve splining diamonds 

and the symmetric line; Ci by the phase boundary, a thin curve splining squares and the symmetric 

line. Schematic spectra of points a-g are shown in Fig. 2. 

FIG. 2. Schematic spectra for the points a-g of Fig. 1. Lines, shaded rectangles and numbers 

denote levels, the continuum and St of the levels or the continuum, respectively. 

FIG. 3. Ground state energy per rung eo vs J± at J x = 0,0.2,0.4,0.6,0.8, and 1. In inset, it vs 

J x at Jj_ = 1. 

FIG. 4. The singlet density per rung p s in the ground state vs J± at J x = 0, 0.2, 0.4, 0.6, and 0.8. 
FIG. 5. Low-lying energies relative to the ground state energy vs J± and J x for the open 

boundary condition. In the Haldane phase, solid triangles linked by solid lines indicate the bottom 

state of the continuum (St = 2). Squares linked by dashed lines the lowest midgap states (MS) 

(St = 1). In the singlet phase, triangles or dots linked by solid lines and diamonds linked by dashed 

lines denote the two lowest excitations with St = 1,0, respectively. Inset for J x = 1. 

FIG. 6. Spin gap vs J x at Jj_ = J» 
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Wang, Figure 2: schematic spectrums 
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Wang, Figure 3: ground state energy per rung 




Wang, Figure 4: density of singlets per rung 




Wang, Figure 5: low-lying excitations 
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